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Abstract 

The spectrum of finite-difference two-magnon operator is investigated for quantum S=l/2 
chain with variable range exchange of the form h(j — k)(x sinh~^ a{j — k). It is found that 
usual bound state appears for some values of the total pseudomomentum of two magnons 
as for the Heisenberg chain with nearest-neighbor spin interaction. Besides this state, a 
new type of bound state with oscillating wave function appears at larger values of the 
total pseudomomentum. 
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The spectral theory of finite-difference operators associated with multimagnon states 
of one-dimensional ferromagnetic Heisenberg chains is far from being complete. From 
previous works, it is known that magnons not only scatter on each other, but can also 
create bound complexes with lower energy. A canonical example of this phenomenon was 
described first by Bethe [1] for the famous case of nearest-neighbor exchange. It was found 
that in this case just one bound state appears for all values of the total pseudomomentum 
P in the interval < P < vr. Much later, it was shown that for next-nearest-ncighbor 
interaction the bound state is always present within this interval but can be destroyed at 
P = TT by second-neighbor interaction [2-4] . The situation for the third neighbor or more 
general nonlocal spin interaction has not been investigated till now. 

In this letter, we shall investigate the two-magnon problem for the integrable ferro- 
magnetic Heisenberg chain with variable range exchange defined by the Hamiltonian 

E hu-k)!^^, (1) 

j^k;j,k=-oo 

where {aj} are Pauli matrices attached to the site j, J > and 

sinh^ a{j — k) ^ 

The range of variation of the parameter a is < a < oo. It will be shown that for all values 
of a in this interval there are two types of two-magnon bound states, the first is similar to 
the bound state for the Heisenberg chain with nearest-neighbor interaction which might 
be considered as a limit of (1-2) as a — > oo (with J oc sinh^ a), and the second corresponds 
to the asymptotically decayed two-magnon wave functions with superimposed oscillations. 

We shall consider two-magnon excitations over the ground state of the model in the 
form 



1 



\^ >= Y: ^{ni,n2)S-^S-^\0 >, S- = -« 



nJi 



where |0 > is the ferromagnetic ground state with all spins aligned up and ip{ni,n2) — 
ip{n2,ni). Acting with (1) on the state \ip >, the equation Hlip >= e\ip > results in 

H^^^ij{ni,n2) = ei>{ni,n2), (3) 
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where the action of the finite-difference operator H^'^^ is given by 

//(2VK,ri2) = -[ Yl {h{ni-k)'ilj{k,n2) + h{n2-k)'ilj{ni,k)) 

+ (2/i(ni-n2)-2£o)V'(ni,n2)], (4) 

where Sq — Y,k^oh{k). A two-parametric solution to the eigenequation (3) with the 
exchange strength (2) reads as follows [5] , 

^(ni,n2) = e^(*'i"i+^'2"^)(coth7(pi,p2) +cotha(ni -712)) 

+e^(f2"i+fi'*2)(coth7(pi,p2) - cotha(ni - n2)), (5) 
where the phase factor 7(]?i,]?2) is expressed via the magnon pseudomomenta pi,2, 

coth7(pi,p2) = 5 ■ (6) 

The formula for two-magnon energy reads 

, , , , , , ^ ^ 2 - cos(n]9i) - cos(np2) 

riT^O ^^^^ l^^J 

where the sum can be expressed via the Weierstrass functions. The structure of the 
function f{p) will be crucial for our treatment. It is given by the expression 

where Ci is the Weierstrass zeta function with quasiperiods 1 and uj — m/a. It admits 
also the representation [6] 



f{p) — cot ^ — a ^ coth -\- an^ + coth — an^ 



(9) 

n=l 

Prom these representations, one can easily find that / is odd and double quasiperiodic in 
P, 

fip) = -fi-p), (10) 

/(p + 27r) = /(p), (11) 
/(p + 2m) = /(p) + 2a. (12) 



The derivative f'{p) is double periodic and due to (11), (12) determined by its restriction 
to the fundamental domain 



D : < 3fJep < 27r, -a < Qmp < a. 



The representations of f'{p) are 



la 



la . ^ 



(13) 



2sin'f 2 sinh^ (I + an) ' 
where pi is the Weierstrass p-function with periods (l,c<;). 

If P11P2 are real, the wave function (5) describes scattering of magnons. In the bound 
state, the wave function must vanish as |ni — ^2] ^00. It means that pi and p2 should 
be complex with P — Pi + P2 real. The simplest possibility is given by the choice 



Pi 



-iq, 



(14) 
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where q is real. One can always choose g > for convenience. Then vanishing of ip{ni, ^2) 
as |ni — 712] — > 00 is equivalent to the condition 



coth7(pi,p2 = 7, = 1- 

Za 



(15) 



Moreover, one can choose q < a. At first sight, the interval a < q < 2a should be also 
taken into consideration but one can see that in this case the rearrangement pi — > p2 + 2ia, 
P2 ^ Pi — 2ia, q ^ 2a — q returns the problem to the interval < g < a due to 
quasiperiodicity property (12). 

The equation (15) under the condition (14) can be rewritten in more detailed form. 



Fp{q) 



1 

2a 



'2tq^ 

TT \2a/ 



I7l\ 



Ci 



iP 



q_ 

2a 



Ci 



iP 



2a 



= 0. 



Aa 2a J \4a 

At fixed real P and real q, the function Fp{q) is real. In the limit a — > 00 it reads 



(16) 



, P iq\ f P iq 

eoti- + | -cot 



e-^-cosP/2 



coshg — cosP/2 

One can see that for each — tt < P < tt there is just one bound state with the solution 
q — — In cos ^ of the equation Fp{q) — 0. In the case of finite a, which is of our interest. 



the situation becomes more involved. One has to investigate zeroes of transcendental 
function Fp{q) (16). To this end, note that (10) and (12) imply the following properties 
ofFp(g), 

Fp(0) = 1, (17) 

Fp{q) = -Fp{2a - g), Fp{q) = -Fp{-q) + 2. (18) 

This means that Fp{a) — 0, i.e. q — a is a, solution to (16) for all P. However, one can 
see by the direct substitution of (14) with q — a into (5) that the wave function in this 
case vanishes identically, i.e. this zero is unphysical and physical solution must lie in the 
interval < g < a. Let us find the existence criterium for it. 

Note first that due to (17), (18) Fp(0) = 1 and Fp{a) = 0. Then a nontrivial zero of 
Fp{q) in the interval < g < a must exist if the derivative of Fp{q) at the point q — a is 
positive, 

'iP 1\ /iP V 



pi 



1\ (ip 

2 j + U 



Aa 2J ° V4a 2 
Since pi{x — 1) = pi (a;), this condition can be cast into the form 

-'/'(f+") = -^<'© + s;«"(s+5) 

1 



> 0. 



~ 2„^^ 



sinh^ ( f + a{n + 1/2) 



> 0. 



(19) 



Let us prove that the inequality (19) has nontrivial solutions. Note first that the 
inequahty definitely holds for P = 0. On the other hand, let us estimate the left-hand 
side at P = TT, 



E 



1 



sinh^ (f + a{n + 1/2 



E 

n=0 



1 



+ 



1 



= -E 



sinh^ (f + a(n + 1 /2)) sinh^ (f - a{n + 1/2 
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< 0. 



^ cosh^(2n + l)a 

It means that there should be some value Pcr(ot), < Pcr{o) < tt, such that inequality 
(19) is satisfied for all values < P < P^, since the left-hand side of (19) is a continuous 
function of P. At these values of the total pseudomomentum there should be at least one 
bound state specified by (14), (15). At P > Per, the inequality (19) does not hold and 



there are no bound states of the type I since the derivative of the left-hand side of (19) 
with respect to P does not change sign in the interval (0, tt) (cf. the known zeroes of 
function p[). 

Let us consider another possibility for getting the bound state. To this end, note that 
Taking into account the relations 



tt" V2ay ' V2. 
one finds that 

/(p + ia) = a + ix{p) , f(p-ia) ^ -a + ix(p) , (20) 
where x(p) is real for real p. Suppose now that 

Pi = Pi + ia, P2 = P2 - ia, (21) 

where both pi and p2 are real, pi ^ p2 (the case pi = p2 corresponds to trivial solution to 
the equation (15) mentioned above). Substitution of (21) into (15) gives with the account 
of (20) 

X(pi)-X(p2) = 0. (22) 



Note also that 



X(0) = X(vr) = 0, 
X'(0) = ? E —^4 A > 0' (23) 



n=— oo 



sinh^ a (n + 



V2J = 2„i:>^sinh^(M + ,(,+ l)) 

It follows from (23) that there should be some value po at which xip) has a maximum 
on the interval [0, tt] and the corresponding Pq = 27r — po at which x{p) has a minimum 
on the interval [tt, 27r]. As a matter of fact, po — Since /'(p) is linearly related to 
the Weierstrass function pi, x'(p) can have only two zeroes in the fundamental domain. 



i.e. there are no other extrema of x{p) the interval [0, tt]. The presence of a maximum 
means that the equation 

xiP) = Xo 

has two distinct real roots if < xo < X {^)^ ^ < Pi < tt and < p2 < These 
roots serve also as nontrivial solution to the equation (22) and thus give, via (21), the 
bound state of type II in which the wave function oscillates and decays exponentially 
as |ni — ^2] — > 00. For Per < P < tt such a solution always exists. Similar solutions 
corresponding to — x < Xo < can be found with any — tt < P < —Per- The 

methods used above do not allow us to say anything definite about the presence of bound 
state at TT < \P\ < 27r different from the just described ones. 

In conclusion, we have proved the existence of two types of bound states of two 
magnons in the model with nonlocal spin interaction (1-2). The states of type I appear 
for small values of the total pseudomomentum P, < \P\ < Per < tt. Due to oscillation 
of terms in the series (7), one can state that the energy for the type I states is lower than 
the continuum boundary only for sufficiently small values of P. The state of the type 
II appears if the total pseudomomentum exceeds the critical value. Per < \P\ < tt. Our 
treatment has been universal with respect to parameter a in the interval < a < 00. In 
the nearest-neighbor limit a — > 00, the type II states with complex relative pseudomo- 
mentum and oscillating wave function disappear {Per ^ tt as it can be seen from (19)) 
and the result coincides with the known one. An interesting open problem concerns the 
presence or absence of multimagnon bound states. The condition of their existence is 
given by the set of transcendental equations of the type (15). At present, we have not 
found the way of their analytical investigation. 
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